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ABSTRACT 

o , 

^SJ ' Aims. In this work we present the study of five 'dying' nearby (z < 0.2) radio galaxies belonging to the WENSS minisurvey and to 

the B2 bright catalogs: WNB 1734-1-6407, WNB1829-f691I, WNB185I-H5707, B20120+33, and B2 1610-^29. 
Methods. These sources have been selected on the basis of their extremely steep broad-band radio spectra, which is a strong indication 
that these objects either belong to the rare class of dying radio galaxies or that we are observing 'fossir radio plasma remaining from 
a previous nuclear activity. We derive the relative duration of the dying phase from the fit of a synchrotron radiative model to the radio 
spectra of the sources. 

Results. The modeling of the integrated spectra and the deep spectral index images obtained with the VLA confirmed that in these 
sources the central engine has ceased to be active for a significant fraction of their lifetime although their extended lobes have not 
yet completely faded away. We found that WNB 1 85 1+5707 is in reality composed by two distinct dying galaxies, which appear 
. , blend together as a single source in the WENSS. In the cases of WNB 1 829+69 II and B2 0120+33, the fossil radio lobes are seen in 

^,1^ ■ conjunction with a currently active core. A very faint core is detected also in a MERLIN image of WNB1851+5707a, one of the two 

dying sources composing WNB 1 85 1+5707. We found that all sources of our sample are located (at least in projection) at the center 
of an X-ray emitting cluster 

Conclusions. Our results suggest that the duration of the dying phase for a radio source in cluster can be significantly higher with 
Q . respect to that of a radio galaxy in the field, although no firm conclusions can be drawn due to the small number statistics involved. 

5— ( ' The simplest interpretation for the tendency of dying galaxies to be found in cluster is that the low-frequency radio emission from 

' the fading radio lobes last longer if their expansion is somewhat reduced or even stopped. Another possibility is that the occurrence 

, of dying sources is higher in galaxy clusters. We argue that radio sources in dense environment, like e.g. the center of cooling core 

clusters, may have a peculiar accretion mode which results in a bursting duty cycle sequence of active and quiescent periods. This 
result could have important implications for theories of the life cycles of radio sources and AGN feedback in clusters of galaxies but 
awaits confirmation from future observations of larger, statistically significant, samples of objects. 

' Key words. Radio continuum: galaxies - Galaxies: active - Galaxies: jets - Galaxies: clusters 

in 

1. Introduction tive losses of the relativistic electrons. The first example of such 

^ ^ . .. , . . .,.„,, sources is B2 0924+30 and was illustrated by Cordey (1987). It 

^ Dymg radio galaxies represent an interesting, but still largely ^j^^ .^^^ ^^^^ ^^^.^ j^^^^^ ^^^^^ intermittently 

O unexplored, stage of the active galactic nuclei evolution. During \u » • » « • i u t » n • » wit 

, , ^ , ' ° ,, , , ° ,,„7 , ° oreventhatets flicker before eventually going off completely. In 

^ their active stage, which may last several 10^ yeai-s, the strong ra- ^^.^ ^^^^^^.^^ ^ ^^^^^^^ ^^^^^j ^^^.^ j^^^^ 

• dio sources associated with elliptical galaxies are supplied with ■ ^^^^ ^^^.^^ ^^^^^ ^j^^^ ^.^^ ^^^^ restarting jets. 

q , energy from active galactic nuclei via plasma beams or jets. Due ^j^^ ^^^^ ^^^^ ^^^^.j ^^^.^ j^^^^ ^^^^ ^.^^ ^ ^^^^.^^^^ ^^^^^ 

;X| . to the continuous accumulation of new partic es, the total spec- 333 ^^^^ ^^^^^^ 2007 for a recent work). The 

^ ; tra of the active radio sources are usually well approximated by ^ ^^^^^^ j^^^^ ^j^^^ ^^^^^^ ^j^^^j disconnected 

, a power law over a wide range of frequencies. The injection of ^^^^ currently active jets. 

energy also sustains the growth of these radio sources which is t .u- 1 i f j- j • .if 

,,,,,, , . , . , In this work we classify a radio source as dying not only if 

governed by the balance between the internal pressure in the f 1 1 u j . u j f a/-'m j\u ■ • 

• ■ ■ • • , the fossil lobes are detached from the AGN and there is no evi- 

radio lobes and the pressure in the hot X-ray emitting external , ^ , 1 \ u . 1 f 1 ■ j f 

, . , , wo , dence for nuclear activity (trivial case), but also if some kind of 

medium in which they must expand (Scheuer 1974). , • . u . f 1 1 u . n j • . 

,■' f .... , . nuclear activity is present but the fossil lobes still dominate the 

At some point, however, the activity in the nuclei stops or , , , ^ j- 1 • m n • u . f n n 

I , ,, I I I n I total source s radio luminosityli Hence, in what follows we will 
falls to such a low level that the plasma outflow can no longer 

be sustained and the radio source is expected to undergo a pe- 1 This choice is motivated by the following two considerations. First, 

nod of fading (dymg phase) before it disappears completely. In ^^^^ tj^e nuclear component we ai'e observing is really produced by 

the dying phase, radio core, well-defined jets and compact hot- new couple of restarting radio jets, we would never know a priori if 

spots will disappear because they are the structures produced by this new activity would be able to last enough to regenerate a radio 

continuing activity. On the other hand, the radio lobes may still galaxy similar to previous event of radio activity whose the fossil lobes 

remain detectable for a long time if they ai-e subject only to radia- represent the remains. Second, it could not be obvious to recognize if 

a nuclear component can be considered as the begin of a new phase of 

Send offprint requests to: M. Murgia, matteo@oa-cagliari.inaf.it AGN activity or rather a fading radio core which is going-off definitely. 
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Table 1. Basic properties for the dying radio galaxies candidates. Spectral indexes and flux densities are from the WENSS and 
NVSS. The spectral curvatures are calculated by measuring aiow between 151 and 325 MHz, for the WENSS minisurvey sources, 
and between 151 and 408 MHz, for the B2 sources. The high-frequency spectral index ahigh is taken between 1400 and 4850 MHz 
for all sources. 



Source 


R.A. (J2000) 


DEC. (J2000) 


redshift 


scale 


S 325MHz 


325MHz 
" 1.4GHz 


SPC 




(h m s) 


(dm s) 




(kpc/") 


(mJy) 




WNB 1734+6407 


17 35 04.6 


+64 06 07.7 


0.1406 


2.4 


120 


2.0 


1.1 


WNB 1829+6911 


18 29 05.6 


+69 14 06.0 


0.204 


3.3 


209 


1.9 


-0.4 


WNB1851+5707a 


18 52 08.5 


+57 1142.3 


0.1068 


1.9 


368 


1.7 


1.1 


WNB1851+5707b 


18 52 09.7 


+57 11 56.6 


0.1066 










B2 0120+33 


01 23 39.9 


+33 15 22.1 


0.0164 


0.33 


727 


1.3 


1.5 


B2 1610+29 


16 12 35.5 


+29 29 05.3 


0.0318 


0.63 


412 


0.9 


0.8 



refer indifferently to dying and re-starting sources if the fossil 
lobes are the dominant component while we shall eventually re- 
serve the term renewed source to indicate the case in which the 
new and the fossil components are roughly comparable in lumi- 
nosity, Uke e.g. in the case of the double-double radio sources 
(Schoenmakers et al. 2000; Saripalli et al. 2002, 2003). 

Indeed, given the comparatively short duration of the radio 
galaxy phenomenon, we could expect a large number of dying 
radio sources. However, only a handful of dying radio galaxies 
in this evolutionary stage are known. According to Giovannini 
et al. (1988) only few percent of the radio sources in the B2 and 
3C samples have the characteristics of a dying radio galaxy. It 
is important to note that these sources represent the last phase 
in the 'hfe' of a radio galaxy. Therefore they must be well dis- 
tinguished from the cluster radio halo/rehc phenomenon, that is 
usually not associated with an individual galaxy (see Ferrari et 
al. 2008 and references therein). 

A possible explanation for the rarity of the dying radio galax- 
ies may be the relatively fast spectral evolution they undergo 
during the fading phase. Synchrotron losses and the inverse 
Compton scattering of the Cosmic Microwave Background pho- 
tons preferentially deplete the high-energy electrons. The fad- 
ing lobes are expected to have very steep (a > 1.3; 5 v oc v"") 
and convex radio spectra characteristic of a population of elec- 
trons which have radiated away much of their original energy 
(Komissarov & Gubanov 1994). In fact, in the absence of fresh 
particle injection, the high-frequency radio spectrum develops 
an exponential cutoff. At this point, the adiabatic expansion of 
the radio lobes will concur to shift this spectral break to lower 
frequencies and the source will disappear quickly. On the other 
hand, if the source expansion is somehow reduced, or even 
stopped, there is still the chance to detect the fossil radio lobe, at 
least at low frequency. 

For the reasons mentioned above, low-frequency selected 
samples such us the Westerbork Northern Sky Survey (WENSS; 
RengeUnk et al. 1997) at 325 MHz and the B2 survey at 
408 MHz (Colla et al. 1970, 1972, 1973) are particularly well- 
suited to search for these elusive fossil radio sources. Parma 
et al. (2007), by cross-correlating the WENSS with the NVSS 
(Condon et al. 1998), discovered six new dying sources and three 
new restarted sources. Other two dying galaxies (the central ra- 
dio source in Abell 2622 and MKW03s) and one possibly restart- 
ing source (MKW07) have been recently found by Giacintucci et 
al. (2007) in a low-frequency survey of nearby cluster of galaxies 
performed with Giant Metrewave Radio Telescope. 

In this work we report on the result of a search for dying 
sources in a WENSS sub-sample (the minisurvey; de Ruiter et 
al. 1998) and in the B2 bright sample described by Colla et 



al. (1975). Both samples have been constructed from the par- 
ent surveys by identifying the radio sources with nearby bright 
galaxies. From these samples we have selected five dying ra- 
dio galaxies candidates whose extremely steep spectra are char- 
acterized by a quasi-exponential drop above a frequency of 
about 1 GHz. These are WNB1734+6407, WNB1829+6911, 
WNB1851+5707, B2 0120+33, and B2 1610+29. 

In order to determine whether these sources were really dy- 
ing objects or fossil lobes associated with restarting radio galax- 
ies, we studied their radio continuum emission with the Very 
Large Array (VLA) at various frequencies. The new observa- 
tions and the analysis of the data already available in literature 
permitted us to derive the detailed integrated radio spectra and 
spatially resolved spectral index images for all the five candi- 
dates. We also had the chance to observe at sub-kpc resolution 
the sources WNB 1851+5707 with the Multi-Element Radio 
Linked Interferometer Network (MERLIN) interferometer. 

Furthermore, we found that all the five sources are unam- 
biguously associated to a diffuse X-ray emission in the Rosat 
All Sky Survey (RASS). In all cases the association is with a 
known cluster of galaxies. 

In this work we report the results of this extensive campaign 
of radio observations. In Sect. 2 we described the criteria of the 
source selection. The descriptions of the VLA and MERLIN ob- 
servations are reported in Sect. 3 and 4, respectively. In Sect. 5 
we present the analysis of the radio spectra. The X-ray environ- 
ment is presented in Sect. 6. Finally, a summary of the work is 
given in Sect. 7. 

Throughout this paper we assume a ACDM cosmology with 
7/o = 71 km s-iMpc-\ Q.„ = 0.3, and Da = 0.7. 

2. Source selection 

2.1. The spectral curvature of dying and restarting sources 

In the active stage the total spectra of the radio galaxies are usu- 
ally well approximated by a power law over a wide range of 
frequencies. Spectral breaks at high frequencies are also often 
observed. The break frequency is due to the radiative losses of 
the synchrotron electrons and it can be related to the magnetic 
field and age of the source. If the integrated spectrum of the 
radio source is dominated by the emission of the radio lobes, 
which accumulated the electrons injected by the jets, it can be 
modeled by a continuous injection model (CI) in which the low- 
frequency spectral index oriow represents the injection spectral 
index of the youngest electron populations, ainj, while the high- 
frequency spectral index is limited to ahigh ^ ^inj + 0.5 at fre- 
quencies greater than the break frequency (see e.g. Murgia et al. 
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1999). As the source ages, the break moves to low frequency 
and the radio spectrum eventually becomes a steep power law in 
which aiow approaches ffhigh- For dying sources however we ex- 
pect an exponential cut-off of the integrated spectrum above the 
break frequency (e.g. Parma et al. 2007). In this case the high- 
frequency spectral index, ahigh, may became much steeper than 
the low-frequency one, aiow. On the other hand, if the activity 
in the galaxy nucleus restarts ahigh can be flatter that aio„. This 
suggests a possible selection criterion based on the radio spec- 
trum curvature to separate active, dying and restarting sources. 
Following the idea of Sohn et al. (2003) we define the spectral 
curvaturq3, S PC, as 

5 fC = (Thigh - ttlow (1) 

As far as the optically thin portion of the radio spectrum is con- 
sidered, in the framework of the CI model the spectral curvature 
will be very close to SPC ~ 0. In fact, the maximum curvature 
we could expect in the active phase is limited to ahigh - ciow ^ 
0.5. 

On other hand, the integrated spectra of dying sources are 
characterized by an exponential cut-off at high frequencies and 
thus for them extreme spectral curvatures (SPC » 0.5) are 
possible. Finally, for restarting sources the injection of new 
electrons will flatten the high frequency spectrum resulting in 
o^high < aiow, and thus S PC < (Murgia, in preparation). 

Obviously, the proposed spectral curvature classification 
must be always considered to have a statistical meaning and 
used by taking into account the following caveats and warn- 
ings. The SPC to be meaningful requires to be calculated over 
a broad frequency window. Because of the nature of the syn- 
chrotron radiation, the spectral curvature of integrated spectra 
generally occurs over a frequency range of at least two order 
of magnitudes. Hence, if the observed spectral window is too 
small, there is always the possibility to miss-classify as active 
a radio galaxy which is indeed dying. Re-starting sources are 
particularly problematic from this point of view since, as the in- 
jection of new radio plasma proceeds, the spectral curvature of 
the source eventually decreases from S PC > 0.5 to SPC < 0, 
so all values of S PC are in principle possible for them (see Sect. 
2.3 for the case of 3C 338). Moreover, individual radio sources 
could not fit in the above picture if their integrated spectrum is 
dominated by Doppler boosted components, like jets, flat spec- 
trum cores, etc. Other common problems in radio surveys, like 
e.g. the blending of unrelated sources or the difficulties in imag- 
ing extended sources, could lead to spurious S PC. For all these 
reasons, the final confirmation of a dying source necessarily re- 
quires the follow-up with deep, spatially resolved, spectral index 
radio imaging of the dying candidate. 

With this strategy in mind (SPC selection and deep spec- 
tral index imaging) we searched for dying and restarting sources 
candidates in two samples extracted from the WENSS and the 
B2 catalogs. 

2.2. The WENSS minisurvey sample 

Dying radio galaxies are more easily detected at low frequency, 
therefore the WENSS at 325 MHz is particularly well-suited to 
search for these elusive objects. We are interested in steep spec- 
trum radio sources identified with bright galaxies as these should 
be prime candidates for dying sources (Parma et al. 2007). The 
WENSS minisurvey contains 402 radio sources identified with 

^ Please note that Sohn et al. 2003 adopted a slightly different defini- 
tion for S PC with respect to the one proposed here. 



either elliptical or spiral galaxies with red magnitude brighter 
than roughly nii- = 17-18. Using cross-correlation with the 
existing catalogs and new observations with the EfFelsberg 100- 
m telescope at 4.8 and 10.4 GHz (Mack et al. in preparation), 
we have obtained spectral information in the frequency range 
from 38 MHz to 10 GHz for about 200 sources brighter than 
5 325 > 30. The spectral curvature can be used to select dying 
radio galaxies candidates in large samples, like the minisurvey, 
where detailed high-resolution images are not available but there 
is a good spectral coverage. 

Three of these 200 sources show a quasi-exponential drop 
of their integrated spectrum with a spectral index as steep 
as 0- > 1.5 at 1.4 GHz. These sources are WNB 1734+6407, 
WNB 1 829+69 1 1 , and WNB 1 85 1 +5707 . 

The spectral curvatures, calculated between 151 and 325 
MHz (for aiow) and between 1400 and 4850 MHz (for ahigh), 
for the three sources are respectively SCP = 1.1,-0.4, and 1.1. 
Thus, WNB 1734+6407 and WNB 185 1+5707 can be classified 
as dying sources candidates, while WNB 1829+691 1 would be a 
restarting source candidate. 

All the three radio sources are associated with elliptical 
galaxies (see Tab. 1 for a list of their basic properties). As we 
will show in the next section, it turn out that WNB 185 1+5707 is 
indeed composed by two distinct dying galaxies that are blended 
together in the WENSS. The two redshifts listed in Tab. 1 have 
been obtained with the DOLORES spectrograph installed at the 
3.5-m Galileo telescope on the Roque de los Muchachos in La 
Palma, Spain (see Parma et al. 2007 for more details on the op- 
tical photometry). 

In order to determine whether these sources were really dy- 
ing objects or relic lobes associated with active radio galaxies, 
we observed their continuum emission in details with the VLA. 
The brighter and more compact of the three, WNB 185 1+5707, 
has been also observed with the MERLIN. 

2.3. The B2 sample 

The B2 bright sample described in Colla et al. (1975) has been 
selected with very similar criteria to the WENSS minisurvey: 
radio sources with a flux density > 200 mJy at 408 MHz have 
been identified with bright nearby galaxies. This complete sam- 
ple counts 54 elliptical galaxies and thus we may expect it may 
contain a few dying sources. For the B2 sample detailed VLA 
images are available. We searched for sources with relaxed ra- 
dio morphologies and extreme spectral curvatures mostly calcu- 
lated between 151 and 408 MHz (for aiow) and between 1400 
and 4850 MHz (for ahigh). Indeed, the prototypical dying source 
B2 0924+30 with S PC ^ 0.8 and the restarting sourc^l 3C338 
(ahas B2 1626+39) with SPC = 0.7 belong to this sample. 

Besides B2 0924+30 and B2 1626+39, we searched for other 
dying or re-starting sources candidates. We found that the radio 
sources B2 0120+33 and B2 1610+29, with a spectral curva- 
ture of respectively S PC -1.5 and SPC = 0.8, can be indeed 
considered dying radio source candidates. The radio source B2 
0120+33 is associated to the efliptical galaxy NGC 507 at the 
center of the Zwicky Cluster 0107.5+3212. The source has a 
weak radio core and the large scale radio emission appears very 
faint and diffuse at a resolution of few arcseconds even at 1.4 
GHz (Parma et al. 1986). The second candidate, B2 1610+29, is 



^ Note that in this restarting source the active core-jets components 
account for just a small fraction of the source's flux density, which is 
dominated by the steep spectrum lobes in the considered frequency win- 
dow. 
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Fig. 1. WNB 1734+6407. Top-left panel image has been taken at 1.4 GHz with the VLA in A+B+C configuration, the restoring 
beam FWHM is 2.9" x 2.8". The contours of the radio intensity are overlaid on the optical DSS2 image. The top-right panel show 
the overlay of the spectral index image between 1.4 and 4.8 GHz with the 4.8 GHz contour levels of the VLA C+D array image. 
The restoring beam FWHM is 5.2" x 5.1". Both the 1.4 and 4.8 GHz images have the same sensitivity level of 10/vJy/beam. Radio 
contours start from a level of 3cr-rms and scale by V2. In the bottom panel we show the integrated spectrum of the source along 
with the best fit of the synchrotron model. 
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Fig. 2. WNB 1829+6911. Top-left panel image has been taken at 1.4GHz with the VLA in A+B+C configuration, the restoring 
beam FWHM is 2.9" x 2.3". The contours of the radio intensity are overlaid on the optical DSS2 image. The top-right panel show 
the overlay of the spectral index image between 1.4 and 4.8 GHz with the 4.8 GHz contour levels of the VLA C-i-D array image. 
The restoring beam FWHM is 5.2" x 4.3". Both the 1.4 and 4.8 GHz images have the same sensitivity level of 13 juJy/beam. Radio 
contours start from a level of 3cr-rms and scale by V2. In the bottom panel we show the integrated spectrum of the source along 
with the best fit of the synchrotron model. 
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Fig.3. WNB 1851+5707. Top-left panel image has been taken at 1.4GHz with the VLA in A+B+C configuration, the restoring 
beam FWHM is 3.7" x 3.3". The contours of the radio intensity are overlaid on the optical DSS2 image. WNB 1851+5707 is 
composed by two distinct dying sources labeled a and h. The top-right panel show the overlay of the spectral index image between 
1.4 and 4.8 GHz with the 4.8 GHz contour levels of the VLA C+D array image. The restoring beam FWHM is 5.3" x 5.0". The 
sensitivity levels of the 1.4 and 4.8 GHz images are 17 and 9yuJy/beam, respectively. Radio contours start from a level of 3cr-rms 
and scale by V2. In the bottom panel we show the integrated spectrum of the source along with the best fit of the synchrotron model. 
The radio spectrum is representative of WNB 185 1+5707 a, the brighter of the two dying sources (see text). 
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u (MHz) 

Fig. 4. B2 0120+33. Top-left panel image has been taken at 1.4 GHz with the VLA in B+C configuration, the restoring beam 
FWHM is 1 1.2" X 10.4". The contours of the radio intensity are overlaid on the optical DSS2 image. The top-right panel show the 
overlay of the spectral index image between 324 MHz and 1 .4 GHz with the 1 .4 GHz contour levels of the VLA C array image. The 
restoring beam FWHM is 16" x 14.5". The sensitivity levels of 324 MHz and 1.4GHz images are 1.4 mJy/beam and 90;uJy/beam, 
respectively. Radio contours start from a level of 3cr-rms and scale by V2. In the bottom panel we show the integrated spectrum of 
the source along with the best fit of the synchrotron model. 
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Fig. 5. B2 1610+29. Top-left panel image has been taken at 1 .4 GHz with the VLA in B+C configuration, the restoring beam FWHM 
is 5.5" X 5.5". The contours of the radio intensity are overlaid on the optical DSS2 image. The top-right panel show the overlay of 
the spectral index image between 1.4 and 4.7 GHz with the 4.7 GHz contour levels of the VLA D array image. The restoring beam 
FWHM is 17" X 17". The sensitivity levels of the 1.4 and 4.7 GHz images are 210 and 80/vJy/beam, respectively. Radio contours 
start from a level of 3cr-rms and scale by V2. In the bottom panel we show the integrated spectrum of the source along with the best 
fit of the synchrotron model. 
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Fig. 6. MERLIN images WNB 1851+5707a at 1.4 GHz. Left, MERLIN+VLA radio contours at 0'.'4 resolution overlaid the DSS2 
image. Contours start at 0.1 mJy/beam and increase by a factor of V2. Right, zoom of the galaxy core. The MERLIN contours at 
0'.'16 resolution are overlaid to the optical DSS2 image. Contours start at 0.012 mJy/beam and increase by a factor of ^J2. WNB 
1851+5707b is not visible in the field of view of the MERLIN observation. 



associated to the nearby cD central galaxy NGC 6086 in Abell 
2162 and has relaxed morphology in the VLA 1.4 GHz image of 
Parma et al. (1986). More recently the source has been imaged 
with the GMRT at 610 MHz by Giacintucci et al. (2007), who 
found a very similar morphology. 

By using all the useful data that are available in the VLA 
archive since then, we studied in detail the spectral properties 
of these objects in order to determine if they really are dying 
sources or not. 



3. VLA observations and data reduction 

3.1. New observations 

We observed the three dying radio galaxy candidates from the 
WENSS minisurvey at 1.4 and 4.8 GHz with VLA in various 
configurations. A summary of this campaign of observations, in- 
cluding the VLA configuration, frequency date and length of ob- 
servations, is reported in Tab. |2] A bandwidth of 50 MHz was 
used for each of the two frequencies. 

Calibration and imaging were performed with the 
Astronomical Image Processing System (AIPS), following 
the standard procedure: Fourier- Transform, Clean and Restore. 
Self-calibration was applied to remove residual phase varia- 
tions. Data from different arrays were combined to improve 
uv-coverage and sensitivity. We combined the A, B and C arrays 
at 1.4 GHz and the C and D at 4.8 GHz. Each combined data set 
was self-calibrated. The relevant image parameters are reported 
in Tab. [3] 

In order to estimate the flux density and the spectral index 
of the three sources the primary beam correction was applied to 
all the images by applying the AIPS task PBCOR. For the pur- 
poses of the spectral index imaging, we convolved the 1 .4 GHz 
images to the slightly larger beam of the 4.8 GHz ones. However, 
we emphasize that the angular resolution of the final images at 
the two frequencies were already well matched, not just by fil- 
tering the data, but by having comparable intrinsic coverage over 



the relevant spatial frequencies. Furthermore, both data sets have 
nearly the same sensitivity levels. 

Table 2. VLA observations summary. 



Source 


Array 


Frequency 


Date 


Duration 






(GHz) 




(hours) 


WNB 1734+6407 


A 


1.4 


21-Apr-2002 


5.2 




B 


1.4 


17-Jul-2002 


2.7 




C 


1.4 


20-Dec-2002 


2.4 




C 


4.8 


21-Dec-2002 


6.2 




D 


4.8 


07-Dec-2001 


3.0 


WNB 1829+69 11 


A 


1.4 


ll-Apr-2002 


2.6 




B 


1.4 


17-Jul-2002 


2.7 




C 


1.4 


20-Dec-2002 


2.4 




C 


4.8 


21-Dec-2002 


0.4 




D 


4.8 


07-Dec-2001 


3.1 


WNB 185 1+5707 


A 


1.4 


26-Mai--2002 


3.0 




B 


1.4 


17-Jul-2002 


0.8 




C 


1.4 


22-Jun-2002 


2.8 




C 


4.8 


15-Dec-2002 


4.0 




D 


4.8 


08-Dec-2001 


3.2 



The VLA total intensity iso-contours of WNB 1734+6407, 
WNB 1829+6911, and WNB 1851+570 at 1.4 and 4.8 GHz are 
respectively shown in top left and right panels of Figs. 1 to 3. 
The 1.4 GHz contours are overlaid to the optical image of the 
red Digital Sky Survey (DSS2), while the 4.8 GHz contours are 
overlaid to the spectral index image between the two radio fre- 
quencies. 

3.2. Archival data 

We recover from the VLA archive all the useful data that are 
available for the radio sources B2 0120+33 and B2 1610+29. A 
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summary of relevant parameters of the images produced is listed 
in Tab. 3. 

The source B2 0120+33 has a useful 325 MHz B-array ob- 
servation in the VLA archive (project AC0845). The observation 
was acquired in hne mode with 15 spectral channels in a band- 
width of 6 MHz per each of the two IFs. The source 3C 48 was 
used as amplitude, phase, and bandpass calibrator. In the final 
imaging the data were averaged to 6 channels and mapped using 
a wide-field imaging technique, which corrects for distortions 
in the image caused by the non-coplanarity of the VLA over a 
wide field of view. A set of overlapping fields was used to cover 
an area of about 2.5° x 2.5° around the radio source. All these 
fields were included in CLEAN and used for several loops of 
phase self-calibration. The central frequency of the final images 
is 324 MHz. At 1.4 GHz we produced a B+C image at a resolu- 
tion of about 10". This image is shown in top-left panel of Fig. 4. 
The resolution of the C array 1.4 GHz image alone well matches 
that of the 325 MHz image and the two have been used to derive 
the spectral index image shown in top-right panel of Fig. 4. It is 
worthwhile to mention that flux densities of both the 325 MHz 
and the 1 .4 GHz images are perfectly consistent with that of the 
WENSS and the NVSS. At 4.8 GHz we were able to combine 
the B and C array. Only the core and the brightest part of the 
lobes are detected in this image that can be used however to de- 
rive an integrated flux density for the radio source. At 8.6 GHz 
we recover a D array image in which only the core of the radio 
galaxy is visible. 

The source B2 1610-H29 has an angular size of about two ar- 
cmin and has a very relaxed morphology which is completely 
resolved out by the A configuration of the VLA at 1.4 GHz. 
Thus, we opted to combine just the B and C array at this fre- 
quency, obtaining a final resolution of 5'.'5 x 5'.'5 and a rms noise 
level of 70 //Jy/beam. This is the highest resolution image we 
present for this source and it is shown as iso-contours in top-left 
panel of Fig. 5. The flux density of our B+C image is consistent, 
within the errors, with the NVSS flux. Therefore, we are confi- 
dent to have recovered most of flux density from the extended 
source structure. We also produced images at 1 .4 and 4.7 GHz in 
C and D configuration, respectively. These images, which have 
been restored with the same resolution of 17" x 17" and cor- 
rected from the primary beam attenuation, are characterized by 
roughly the same uv-coverage and have been used to produce 
the spectral index image shown in the top-right panel of Fig. 5. 
Finally, for the purposes of the analysis of the source integrated 
spectrum, we reduced separately the two IFs of the L band obser- 
vation obtaining a measure of the source flux densities at 1452 
and 1502 MHz. 

3.3. Comments on individual sources 

The VLA integrated flux density and spectral index for the indi- 
vidual source components are hsted in Tab. 4. 

3.3.1. WNB 1734+6407 

Characterized by two relaxed lobes lacking hot-spots, the radio 
morphology of WNB 1 734+6407 resembles that of B2 0924+30, 
which is considered the prototype of fossil radio galaxies 
(Cordey 1987, but see also Jamrozy et al. 2004 for more recent 
data on this source). The radio spectrum of the fossil lobes of 
WNB 1734+6407 is so steep that their surface brightness at 6 
cm is below the sensitivity level of our observations. Thus, we 
can only place lower hmits on the spectral index which result 



Table 3. Image parameters summary. 



Source 


Array 


V 


Beam 








(GHz) 


(arcsec) 


(//Jy/beam) 


WNB 1734+6407 


A+B+C 


1.4 


2.9 x 2.8 


10 




c+u 


A O 

4.6 


5.Z X 5.1 


1 A 

10 


WiNDloZb'+Db'l i 




1 A 
1.4 


L.y X Z.J 


1 'i 




C+D 


4.8 


5.2 X 4.3 


13 


WNB 185 1+5707 


A+B+C 


1.4 


3.7 X 3.3 


17 




C+D 


4.8 


5.3x5.0 


9 


B2 0120+33 


B 


0.324 


16.0 X 14.5 


1400 




C 


1.4 


13.4 X 12.5 


50 




B+C 


1.4 


11.2 X 10.4 


80 




B+C 


4.8 


4.0 X 4.0 


30 




D 


8.4 


8.3x5.1 


70 


B2 1610+296 


B+C 


1.4 


5.5x5.5 


70 




C 


1.4 


17.0 X 17.0 


210 




D 


4.7 


17.0 X 17.0 


80 



in a > 2.7 and a > 2.3 for the northern and southern lobe, re- 
spectively. We also observe a slightly extended component coin- 
cident with the galaxy center. Also this feature has a quite steep 
spectral index, a = 1.3, and its nature remains unclear. In the 
VLA A-array image (not shown) this feature is elongated in N-S 
direction with an extent of about 7 kpc x 14 kpc, resembling a 
core-jet morphology. 

3.3.2. WNB 1829+6911 

The radio appearance of WNB 1829+69 11 is virtually identi- 
cal to that of 3C 338, a nearby radio source associated with 
the central dominant galaxy in the coohng flow cluster AbeU 
2199 (Jones & Preston 2001). In both sources we observe the 
presence of fossil plasma remaining from a previous activity 
in conjunction with a restarting core. The extended emission of 
WNB 1 829+69 1 1 has a very steep spectrum, or > 1.6 and a > 1.9 
respectively for the W and E-lobe. The source exhibits a bright 
core with a much flatter spectrum, a = 0.6. Most probably this 
region of the radio source is powered by a new couple of restart- 
ing jets. The core emission is responsible for the high-frequency 
flattening seen in the source integrated spectrum (Fig. 2 bottom 
panel). An arc-like feature with a very steep spectrum (a > 2.5) 
can be seen on the SW side of the core (Fig. 2 top-left panel). 

3.3.3. WNB 1851+5707 

Our VLA images reveal that WNB 185 1+5707 is in reality 
composed by two distinct dying radio galaxies. For source 
WNB1851+5707a we measured a spectral index a = 2.2 while 
WNB1851+5707b we measure a - 2.3. This association is a re- 
ally intriguing fact considering the rarity of this kind of sources. 
According to the interpretation that these are relic sources, we 
found that, at the sensitivity hmits of our VLA observations, they 
possess neither radio jets on kpc scale nor bright hot spots inside 
the fossil lobes. The two hosting galaxy have nearly the same 
redshifts (see Tab. 1), making the possibility of a spatial casual 
coincidence even more unlike. 

It should be noted fliat WNB 1 85 1 +5707 is located at the cen- 
ter of a galaxy cluster, see Sect. 6. Thus, it could be that the two 
dying lobes of WNB1851+5707b are indeed two previous out- 
bursts of WNB1851+5707a that are emerging from the cluster 
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Distance along the lobes (kpc) 

Fig. 7. B2 0120+33 spectral index profile. Each point represents the average over a circular region as shown in the inset. The solid 
line is the best fit of the C/qff rnodel described in the text. The reference dashed lines correspond to different relative durations of 
the dying phase. 



center because of the buoyancy forces. Although this scenario 3.3.5. B2 1610+29 
seems less likely, in view of the almost perfect spatial coinci- 
dence of WNB1851+5707b with the second galaxy. 



3.3.4. B2 0120+33 



The radio source is asymmetric, the west lobe being brighter and 
smaller that the east one. The radio emission is very weak and 
diffuse. A similar morphology has been observed in GMRT im- 
ages at 235 and 610 MHz by Giacintucci et al. 2010 (in prepara- 
tion). At arcsecond resolution the only compact structure is the 
faint unresolved core. There is no evidence of kpc-scale jets in 
any of our images. Hence, the extended lobes of B2 0120+33 
have been produced by an earlier phase of activity. 



The morphology of B2 1610+29 is exactly the one we would 
expect for a dying radio source: two symmetric and relaxed ra- 
dio lobes sit on the opposite side of the host galaxy, there are 
no jets nor hot-spots in the lobes, and the host galaxy lacks a 
bright radio core. All these features are strong indications that 
B2 1610+29 is a dying source where only the fading lobes are 
still visible. The source has an overall linear size of about 100 
kpc with a small offset of the source major axis to south-est 
with respect to the position of the galaxy. The host galaxy of 
B2 1610+29 has been recently studied with the VLBA at 5 GHz 
by Liuzzo et al. (2010). They did not detected any compact ra- 
dio structure in their high resolution images. Also this results 
suggests that most likely B2 1610+29 is a dying radio source. 
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Fig. 8. B2 1610+29 spectral index profile. Each point represents the average over a rectangular region as shown in the inset. The 
solid line is the best fit of the C/qff rnodel described in the text. The reference dashed lines correspond to different relative durations 
of the dying phase. 



4. MERLIN observations and data reduction 

Among the five dying sources, WNB1851+5707a is the most 
compact and bright. We observed this source with the MERLIN 
interferometer in order to unveil the possible presence of a core 
or jets unresolved by the VLA images. WNB1851+5707a was 
observed in three days on April 2005 for a total on-source time 
of about 18 h with MERLIN (6 antennas) and Lovell. The ob- 
serving frequency was 1.408 GHz, with a bandwidth of 16 MHz 
in both circular polarizations; the data were taken in spectral-line 
mode (32 x 0.5-MHz channels). OQ208 (1.00 Jy) and 3C286 
(14.79 Jy) were used as flux calibration sources, while 1851+609 
(0.25 Jy) as a phase calibrator. The 1.4 GHz data from MERLIN 
and the VLA A-array were combined to increase the brightness 
sensitivity and uv-coverage, in order to image finer details of the 
components seen in the L-Band MERLIN images, whilst also 
imaging the more diffuse radio emission that is observed in the 
VLA images. The length of the shortest baseline of the MERLIN 



array at 1 .4 GHz is 6 km. At this frequency the interferometer 
is then not sensitive to structures larger than fJ'na™^ - (^-^ 
kpc at the distance of WNB 185 1+5707 a). The combination of 
the 1.4 GHz MERLIN data with the VLA A-array data gives a 
shortest baseline of 0.68 km, so that ^l^tx^^ - 64" (136 kpc). 
Several images were produced using the AlPS task IMAGR, 
and deconvolved with the multi-scale CLEAN algorithm. The 
full-resolution 0'.'16 MERLIN image is shown in the right panel 
of Fig. 6. The only detected feature is a faint point-like source 
with a flux density of 0.4 ±0.1 mJy at the center of the op- 
tical galaxy. There is no evidence of active radio jets. Indeed, 
most of the flux density of WNB 1851 +5707a comes from the ex- 
tended structure already mapped in the VLA observations, con- 
firming that WNB1851+5707ais a "de-energized" radio source. 
The MERLIN image is presented in combination with the VLA 
A-array at a resolution of 0'.'4 in the left panel of Fig. 6. The 
MERLIN+VLA contours reveal that WNB1851+5707a is an 
amorphous radio source with a linear size of roughly 30 kpc x 
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Source 


Component 


size 
kpcxkpc 


S 1.4 

mJy 


5' 4.8 

mJy 


"l4 


WNB 1734+6407 


N-Lobe 
o-Luue 
Core-Jet? 


63x42 

DU A +Z 

14x7 


2.9 + 0.03 

1 7-1-0 0^ 
i . / It U.Uj 

0.7 ±0.01 


< 0.1 
^ n 1 

< U. 1 

0.14 ±0.02 


> 2.7 

> Z..J 

1.3 ±0.1 




Total 


120 X 62 


5.9 ±0.1 


0.14 ±0.02 


3.0 ±0.1 


WNB 1829+69 11 


E-Lobe 
W-Lobe 
Arc 
Core 


72 X 47 
69x41 
33 X 16 


2.1 ±0.07 

1.4 ±0.06 
2.1 ±0.05 

4.5 ± 0.02 


< 0.2 

< 0.2 
<0.1 

2.1 ±0.03 


> 1.9 

> 1.6 

> 2.5 
0.62 ±0.01 




Total 


135 X 61 


11.2±0.1 


2.1 ±0.03 


1.4 ±0.01 


WNB1851+5707a 


Total 


15 X 30 


41.8 ±0.08 


2.8 ± 0.04 


2.2 ±0.1 


WNB1851+5707b 


N-Lobe 
S-Lobe 


23 X 17 
19 X 19 


2.6 ± 0.07 
2.3 ± 0.05 


0.17 ± 0.02 
0.13 ±0.02 


2.2 ± 0.3 

2.3 ±0.3 




Total 


50x25 


4.9 ± 0.06 


0.3 ± 0.03 


2.3 ±0.2 


B2 0120+33 


E-Lobe 
W-Lobe 
Core 


47 X 24 
32 X 22 


43.6 ±0.3 
56.5 ±0.2 
2.1 ±0.2 


3.64 ± 0.6 
4.5 ±0.5 
1.8±0.1 


2.0 ±0.14 

2.1 ±0.1 
0.13 ±0.09 




Total 


73x40 


99.5 ± 1.5 


8.6± 1.1 


2.0 ±0.1 


B2 1610+29 


Total 


100 X 30 


100 ± 1.4 


15.5 ±0.3 


1.5 ±0.02 



15 kpc. The radio source is also considerably asymmetric being 
much more extended to the south of the optical galaxy. 



5. Analysis of the radio source spectra 

5.1. Compilation of integrated spectra 

In addition to our own observations, we collected all the spectral 
information available in the literature for the five dying radio 
galaxies. We made use of the CATS (the on-line Astrophysical 
CATalogs support System, at http://cats.sao. ru/; Verkhodanov et 
al. 1997) to recover data from catalogs at different frequencies. 
We put all the flux densities from all these reference on the ab- 
solute flux density of Baars et al. (1977) by scaling for the mul- 
tiplicative factor Usted in Helmboldt et al. (2008). The resulting 
catalogued flux densities given in Col. 3 of Tab. 5. 

The various catalogs we used to compile the integrated spectra 
have quite different angular resolution. It is therefore necessary 
to remove the flux density due to confusing sources in order to 
recover a reliable spectrum for each dying source. We estimated 
the contribute of the confusing sources at each frequency on the 
basis of the VLA images at 1.4 and 4.8 GHz. The confusion flux 
density at a given frequency is given by the sum 



^ confusion 



= £^'1425(^1425)- 



(2) 



/■=i 



where v is the frequency in MHz, S 1425 is the flux density mea- 
sured in the VLA image at 1425 MHz, a, is the spectral index 
of the confusing source measured from the VLA images at 1 .4 
and 4.8 GHz, and Ny is the number of sources that fall within the 
beam of the catalog at frequency v. Finally, we calculated the 
source flux density at a given frequency as 



this is the quantity listed in Col. 7 of Tab. 5. 

The total spectra of the five dying radio sources are 
show in bottom panel of Figs. 1 to 5. For four of them, 
namely WNB1734+6407, WNB1851+5707, B2 0120+33, and 
B2 1610+29, the integrated spectrum presents a strong expo- 
nential cutoff in the observed frequency range. However, WNB 
1829+6911 shows an evident flattening of the integrated spec- 
trum at high frequency. This flattening is due to the presence of 
the core-jet component that dominates the spectrum at the high- 
est frequencies. 

5.2. Spectral modeling of the integrated spectra 

We modeled the integrated spectra assuming the radiative en- 
ergy losses to be dominant with respect to other processes (e.g. 
adiabatic losses). The pitch angles of the radiating electrons are 
assumed to be continually isotropized in a time that is shorter 
than the radiative time scale. According to this assumption the 
synchrotron energy losses are statistically the same for all elec- 
trons. After its birth the source is supposed to be fuelled at a 
constant rate (i.e. the continuous injection phase) by the nu- 
clear activity, for a duration fci- The injected particles are as- 
sumed to have a power law energy spectrum A^(e) oc e"*"j, which 
will result in a power law radiation spectrum with spectral in- 
dex Qfinj = (dinj - l)/2. In this phase the source radio spectrum 
changes as a function of time in a way described by the shift of 
break frequency Vb to ever lower values as the time, f „ increases: 



B 



(B2 



(4) 



c _ r- catalogued r. confusion 

O V — O y O y 



(3) 



where B and B/c - 3.25(1 + z)^ are the source magnetic field 
and the inverse Compton equivalent magnetic field, respectively. 
Below and above Vb the spectral indices are respectively ainj and 

Q'ini+0.5. 
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At the time fci the power supply from the nucleus is 
switched-ofF. After that a new phase of duration foFF begins (i.e. 
the dying phase). A new break frequency Vhhigh then appears, 
beyond which the radiation spectrum drops exponentially. This 
second high frequency break is related to the first by: 



ythigh - Vb 



(- 



(5) 



where ts = fci + ^off is the total source age (see e.g. 
Komissarov & Gubanov 1994, Slee et al. 2001, Parma et al. 
2007). 

Thus, the above synchrotron model (hereafter CIqff) is de- 
scribed by four parameters: 

i) ainj, the injection spectral index; 

ii) Vb, the lowest break frequency; 

iii) toFFlh, the dying to total source age ratio; 

iv) norm, the flux normalization. 

In the CIoFF model the magnetic field strength is assumed to 
be uniform within the source. 

The fit of the CIqff is shown as as line in bottom panels of 
Figs. 1 to 5 while the best fit parameters are listed in Tab 6. The 
fits are very good for all the five dying sources. 

- In the case of WNB 1734+6407 the overall radio spectrum is 
particularly steep with an exponential cut-off' beyond a fre- 
quency of about 230 MHz. The spectral fit indicates that the 
source spent in the dying phase about 30% of its total age 
(fOFF/fs = 0.3). 

- The source WNB 1829+691 1 is a different case. Its low fre- 
quency spectrum is rather steep, but strongly flattens above 
1 .4 GHz, due to the presence of a bright core that is clearly 
detected in the 4.8 GHz images. The core have a flat spec- 
trum characterized by a = 0.6 between these two frequencies 
and is likely still active. In order to account for this compo- 
nent in the spectral fit, we added a power law to the CIqff 
model, using the observed spectral index and flux density as 
the normalization. In this case we may be observing fading 
lobes (produced by a previous duty cycle), in conjunction 
with restarting activity in the core. However, we cannot con- 
straint from the spectral fit the relative duration of the dying 
phase since the rising of the new flat spectrum core canceled 
out the the second high break frequency Vhhigh of the fossil 
lobes. 

- Both dying sources a and b contribute to the total spectrum 
of WNB 1851+5707. However, WNB 1851+5707a is ten 
times brighter of WNB 1851+5707b (see Tab. 4). Indeed, 
if the radio sources have a similarly curved spectral shape, 
we can reasonably assume that the total spectrum of WNB 
1851+5707 is representative of source a alone. The radio 
spectrum of WNB 1851+5707a is consistent with an injec- 
tion spectral index of about o-inj = 0.5 and cut-off almost ex- 
ponentially beyond a break frequency of 220 MHz. The du- 
ration of the dying phase is 20% of the total source age. For 
what concerns WNB 1851+5707b we do not have the possi- 
bility to measure its complete total radio spectrum, given its 
closeness to its brighter companion. We can just report the 
particularly steep spectral index of the two fading lobes as 
seen in our VLA images: ffj-^ ^ 2.3. 

- We fitted the integrated spectrum of B2 0120+33 with the 
CIqff model plus a power law which accounts for the emis- 
sion of the radio core. The radio core however is very 



faint and its emission becomes important only at frequen- 
cies greater than 10 GHz. At lower frequencies the source 
spectrum is dominated by the fading radio lobes and it can 
it be modeled with an injection spectral index of Cinj = 0.6 
and a break frequency of about 440 MHz. According to the 
spectral fit the duration of the dying phase should be about 
30% the source total age. 
- Finally, the case of B2 1610+29. The source has an injection 
spectral index of ainj = 0.6 and a break frequency some- 
what higher with respect to the other dying source presented 
above, v/, ^ 5200 MHz. Due to the lack of spectral infor- 
mation at frequencies higher that 10 GHz, we can only put a 
lower limit on the duration of the dying phase: fQFF/fs > 0.4. 
However, for this source we have a very good spatially re- 
solved spectral index image and we can try to get some addi- 
tional constraints on the source age from the fit of the spec- 
tral index profile. 



5.3. Modeling of the spectral index profile in B2 0120+33 and 
B2 1610+29 

The two dying galaxies B2 0120+33 and B2 1610+29 are ex- 
tended enough to permit a study of the spectral variation along 
the fading lobes. 

For B2 0120+33 the profile of the spectral index between 
325 MHz and 1.4 GHz can be traced all along the lobes, see 
Fig. 7. The spectral index trend has been obtained by averag- 
ing the flux densities at the two frequencies in circular regions 
of 36" in diameter centered as shown in the inset of Fig. 7. The 
regions are much larger than the beam so that the spectral index 
measurements are effectively independent. The observed spec- 
tral index trend is rather smooth with a^^ close to about 1 .2 over 
the whole radio source. Assuming a constant source expansion 
velocity the spectral index at a given distance from the core can 
be related to the synchrotron age of the electrons at that location. 
We assume that in B2 0120+33 the youngest electrons were in- 
jected by the jets close to the core during the active phase. This 
behavior is the typical one we observe in many FRI-type tailed 
radio sources (see e.g. Parma et al. 1999). In this case, we expect 
that the break frequency scales as: 



vt{d) = 



Vi,(<imax) 



[(-i/t/max) ■ (1 - toFF/ts) + tOFF/ts? 



(6) 



where the distance d ranges from d - at the location of 
the galaxy, up to - d^ax at the edge of each lobe. Thus, the 
break frequency along the lobes varies from a minimum of vt,, at 
d - c/max, up to a maximum of Vbhieh, at c/ = 0. The two limiting 
break frequencies Vb and Vbhigh are exactly the same as given in 
Eqs. ^ and (|5]l, respectively. In Fig. 7 we report the fit to the 
observed trend obtained in the case of o-inj = 0.65 for different 
values of foFF/fs- It is clear that low values of Tqff/^s correspond 
to spectral index gradients which are too steep if compared to 
the observed one. As foFF/^s increases, however, the expected 
spectral index trends get flatter and flatter and a best fit is found 
for fQFF/fs > 0.8. 

In the case of B2 1610+29 the profile of flie 1.4 to 4.8 GHz 
spectral index can be traced all along the lobes, see Fig. 8. The 
spectral index trend has been obtained by averaging the flux den- 
sities at the two frequencies in boxes perpendicular to the source 
major axis. The boxes are wide as the beam so that the spec- 
tral index measurements are effectively independent. The ob- 
served spectral index trend is rather smooth with ranging 
from about 1.2 at the edge of the lobes up to about 1.4 at the 
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Fig. 9. RASS X-ray images with overlaid the VLA C-array contours at 1 .4 GHz. From top-left to bottom-right: Abell 2276, ZwCl 
1829.3-H6912, RXC J1852. l-)-571 1, ZwCl 0107.5+3212, and Abell 2162. Radio contours stait from a 3o- rms level and increase by 
a factor of V2. 
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center of the radio source. In this particular case, we assume that 
the electrons were injected by the jets close to the end of each 
of the two lobes during the active phase. This spectral behaviour 
is common in FRII type radio sources but is also found in some 
FRI type radio sources (see Parma et al. 1999). In this case, we 
expect that the break frequency scales as: 



Vfe(<imax) 



{[(dmm. - d)/dmaK] ■ (1 - toFp/ts) + tow/ts}^ 



(7) 



where the distance d ranges from d = at the location of 
the galaxy, up to J = <iniax at the edge of each lobe. Thus, the 
break frequency along the lobes varies from a minimum of Vb, at 
d = 0, down to a maximum of Vbhigh, at J = Jmax- In Fig. 8 we 
report the fit to the observed trend obtained in the case of ainj - 
0.6 for different values of topp/ts- It is clear that low values of 
toFp/h correspond to spectral index gradients which are too steep 
if compared to the observed one. As foFp/fs increases, however, 
the expected spectral index trends get flatter and flatter and a best 
fit is found for toHs/h = 0.8 ± 0.1. 

The behavior of the spectral index seen in B2 0120+33 and 
B2 1610+29 could be typical of many dying sources. In practice, 
as a result of the jets switch-off, any pre-existent spectral index 
gradient along the lobes is quickly canceled since the break fre- 
quency reaches roughly the same value in each part of the source. 
We can therefore expect that extreme dying sources, those for 
which foFp/fs - 1, are characterized by very uniform spectral 
index distributions along the fading lobes. As the dying source 
gets older, the spectral index increases systematically but with 
small variations from point to point. 

5.4. Equipartition parameters 

We performed the minimum energy calculation for the five dying 
galaxies by considering the source power at 15 1 MHz, where the 
energy losses of the synchrotron electrons are less dramatic, and 
the volumes as measured from the VLA images at arcsecond 
resolution. In the calculation of the sources power we made use 
of the 151 MHz flux density after the deconvolution of confusing 
sources. The resulting monochromatic powers in the sources rest 
frame, L151, are listed in Tab. 6 and place them below the FRI- 
FRII division (Fanaroff & Riley 1974). 

We assume that the radio sources contain relativistic par- 
ticles and magnetic fields uniformly distributed and in energy 
equipartition conditions. The equipartition parameters (magnetic 
field B mm and energy density Wmin) are generally computed as- 
suming that the relativistic particle energies are confined be- 
tween a minimum eio„ and a maximum ehigh, corresponding to 
the observable radio frequency range, typically 10 MHz - 100 
GHz (see, e.g., Pacholczyk 1970). This choice minimizes the 
source energetics required by the observed radiation in the radio 
band. However, a fixed frequency range corresponds to an en- 
ergy range that depends on the source magnetic field, which may 
change from source to source. A fixed frequency range com- 
putation would miss the contribution from lower energy elec- 
trons, since the e/ow corresponding to 10 MHz is larger than 200 
MeV for B^i^ < 30 fiG (Brunetti et al. 1997 and Beck & Krause 
2005). Because of this we have computed the equipartition pa- 
rameters assuming a fixed low energy cutoff eiow = 10 MeV. The 
high energy ejugh cut-off is chosen to match the spectral break 
in the emission spectrum, Vb, since no electrons are present in 
the source beyond this hmit, see Parma et al. (2007) for further 
details. 



The equipartition parameters computed in this way are re- 
ported in Tab. 6. A comparison with fixed frequency range cal- 
culations shows that the values of our B,^ are larger by up to a 
factor two. 



5.5. Radiative ages 

Assuming a constant magnetic field and neglecting expansion 
losses, the total source age can be calculated from the break fre- 
quency, Vb'. 



u = 1590 



B' 



.0.5 



(B2 + B2^)[(l+^)Vb]0-5 



(8) 



where the synchrotron age fj is in Myr, the magnetic field in 
/iG, the break frequency Vb in GHz, while the inverse Compton 
equivalent field (see Sect. 5.2). 

By adopting the equipartition value Bmin for the magnetic 
field strength we can thus derive the synchrotron age fs. Finally, 
from the ratio which is also given by the fit, we can de- 

termine the absolute durations of the active and dying phases, 
fci and foFF- The total ages of the five dying sources are of the 
order of 10** yrs (with uncertainties up to 50%). For all sources 
but WNB 1829+6911 we are also able to estimate the absolute 
duration of the active and dying phases that result in the range 
10^ - 10^ yrs. We note that these numbers should be consid- 
ered us upper limits to the real source ages. It is likely that ex- 
pansion losses may have played a role during the early stages 
of the radio source growth, and, in the same manner, the mag- 
netic field strength was presumably higher than the actual values. 
Considering both these effects, the shift of Vb to low frequencies 
may have proceeded much rapidly than expected by the simple 
radiative model we considered, and hence, the ages would have 
been overestimated. 

6. The X-ray environment 

The gaseous environment in which radio galaxies are embed- 
ded may play a fundamental role in the latest stages of the radio 
source life. Although no firm conclusions can be drawn about the 
environment of dying sources, due to the small number of ob- 
jects involved, it seems that there is a tendency for dying sources 
to reside in dense enviroimients. Parma et al. (2007) found that 
about half of the dying sources of their sample are located in 
clusters of galaxies while only a few appear to be isolated. 

We have made a search in the RASS for possible X-ray coun- 
terparts of our five dying radio galaxies. We extracted the 0.1 - 
2.4 keV image of a 15' X 15' field around each radio source. The 
X-ray images were corrected for the background and smoothed 
with a (T = 45" Gaussian kernel. The RASS count-rate im- 
ages are shown in Figs. 9 along with the 1 .4 GHz VLA C-array 
iso-contours overlaid. The five dying galaxies presented here do 
not represent an exception to the rule: each source is located, 
at least in projection, at the center of an X-ray emitting clus- 
ter or galaxy group. The clusters of galaxies are Abell 2276, 
ZwCl 1829.3+6912, RX J1852.1+5711, ZwCl 0107.5+3212, 
and Abell 2162 for the radio sources WNB 1734+6407, WNB 
1829+691 1, WNB 1851+5707, B2 0120+33, and B2 1610+29, 
respectively. 

6.1. Abell 2276 

Very little is known about this cluster of galaxies besides that 
Jones & Forman (1999) report an X-ray luminosity 1.02 x 10** 



M. Murgia et al.: Dying Radio Galaxies in Clusters 



17 



erg/s in the 0.5-4.5 keV band based on Einstein data (all X-ray 
luminosities have been rescaled to our adopted cosmology). The 
dying source WNB 1734+6407 is associated with the central 
galaxy at the peak of the X-ray emission, see top-left panel in 
Fig. 9. 

6.2. ZwC1 1829.3+6912 

This galaxy cluster is also known as MACS J 1829.0+69 13 or 
RX J1829. 0+6912. Only sparse information are present in liter- 
ature for this cluster too. It appears as catalog number NEP5700 
in the ROSAT North Ecliptic Pole survey by Gioia et al. (2003). 
These authors report an X-ray luminosity of 1.06 x lO"*"* erg/s 
in the 0.5-2.0 keV band. WNB 1829+691 1 associated to the cD 
galaxy located close to the X-ray peak, see top-right panel in 
Fig. 9. 

6.3. RXJ1 852.1 +5711 

The extended ROSAT Bright Cluster Sample (Ebeling et al. 
2000) lists an intrinsic X-ray luminosity in the 0. 1-2.4 keV band 
of 0.95 X 10** erg/s for this cluster The dying sources WNB 
1 85 1 +5707a is associated to the galaxy at the center of the X-ray 
emission while the dying source WNB 1851+5707b and its host 
galaxy he 30"to the north east. WNB 1851+5707a and WNB 
1851+5707a appear blended together at the resolution of the 
VLA C-array shown as iso-contours in the middle-left panel of 
Fig. 9. 

6.4. ZwCI 0107.5+3212 

The dying source B2 0120+33 is identified with the galaxy NGC 
507 at the center of the Zwicky cluster 0107.5+3212 and is one 
of the brightest galaxies in a very dense region. There are 80 
galaxies with magnitude nipg < 15.7 within a projected distance 
of 1.5 Mpc from NGC 507. This is the so-called Pisces Group 
(Parma et al. 1986). The X-ray emission has been studied in de- 
tail by many authors with different satellites, see e.g. Kraft et al. 
(2004) and references therein. NGC 507 is surrounded by an ir- 
regular large-scale X-ray emission with a peak coincident with 
the position of the optical galaxy, and several secondary peaks, 
due to the interaction between the radio lobes and the surround- 
ing gas (Paolillo 2003). Kim & Fabbiano (1995) report an X-ray 
luminosity of 0.8 x lO*^ in the 0.1-2.4 keV band for the NGC 
507 group. Kraft et al. (2004) observed a sharp edge in the X-ray 
surface brightness profile in a Chandra observation of NGC 507. 
These authors suggested that the discontinuity is likely caused 
by an elemental abundance jump driven by the inflation of the 
radio lobes and entrainment of material from the central regions 
of the cluster. 



6.5. Abell2162 

Abell 2162 is member of the Hercules super cluster complex 
(Einasto et al. 2001). The X-ray emission of this nearby cluster 
of galaxies has been studied by Ledlow et al. (2003) who report 
an X-ray luminosity in the 0.5-2.0 keV band of 0. 16 x 10*^ erg/s 
based on ROSAT data. The fading radio lobes of B2 1610+29 
are associated to NGC 6086, the brightest galaxy of the cluster 
located close to the X-ray peak (Fig. 9 bottom panel). 



Indeed, these results suggest that the occurrence of dying 
sources could be higher in galaxy clusters, confirming the find- 
ings of Parma et al. (2007). 

In order to determine if the environment really plays a role 
in increasing the occurence or the duration of the dying phase 
we need to estimate the relative fraction of dying sources in iso- 
lated and cluster galaxies. This is not trivial because estimating 
the completeness of the inhomogeneous radio, optical, and X- 
ray samples is hard. However, by imposing safer selection limits 
we can attempt to perform some very rough statistical analysis. 
A complete sample of 1 19 sources can be constructed from the 
minisurvey by considering all galaxies with optical magnitudes 
nir < 16 and radio flux density 5 325MHz > 30 mJy (de Ruiter 
et al. 1998). If we further restrict our attention to the elliptical 
galaxies we have a complete sample of 90 radio sources. Note 
that we consider WNB 1851+5707 as two distinct radio galax- 
ies, both in the dying phase. 

For each of these radio sources we searched in the NED 
database for the presence of a galaxy cluster or group to within 
10' from the position of host galaxy. We found that 21 out of 
90 sources lies, in projection, in Abell, Zwicky or X-ray se- 
lected galaxy clusters or groups, while we considered the rest 
to be isolated. In this complete sample there are three of our dy- 
ing sources: WNB 1734+6407, WNB 1851+5707a and WNB 
1851+5707b. The restarting source WNB1829+6911 exceeds 
the optical limit due to its larger distance and therefore is not in- 
cluded in the complete sample. Therefore, the fraction of dying 
source in general is ~ 3.3% (3/90). However, this fraction de- 
pends on the environment since the frequency of dying sources 
in clusters is 14% (3/21). 

In the complete B2 bright sample (Colla et al. 1975), we 
found the presence of fossil radio lobes in four of 54 ellipti- 
cal galaxies. Of these, three (B2 0120+33, B2 1610+29, and 
B2 1636+39) are in cluster and only one (B2 0924+30) can be 
considered relatively isolatecQ. With the same criterion describe 
above, we estimated that about 21 radio sources of the B2 sample 
are located in projection close to the center of a galaxy cluster 
or group. The fraction of dying source in general is 7% (4/54), 
while the fraction of dying sources in clusters is 14% (3/21). 

The fractions of dying sources (both in general and in clus- 
ter) in the minisurvey and in the B2 samples are consistent to 
within the large uncertainties. We then merged the statistics of 
the two samples in order to estimate the probability for a dying 
source to be in cluster: 

P{dying\cluster)P{cluster) 

P{cluster\dying) = — — — , (9) 

Pyaymg) 

where P{dying\cluster) - 6/42 is the fraction of dying sources in 
cluster, Picluster) - 42/144 is the total fraction of radio sources 
in cluster, and P(dying) - 7/144 is the total fraction of dying 
sources. Indeed, we obtained that the probability for a dying 
source to be in cluster is P(cluster\dying) - 6/1, i.e. ~ 86%. 
We estimated through a simple Monte Carlo simulation that if 
the possibility to have a dying source were independent from 
the environment, the probability to have 6/7 of dying sources in 
cluster by chance is less than 0.5%. 

One possibility is that the low-frequency radio emission 
from the fading radio lobes last longer if their expansion is some- 
what reduced or even stopped. Another possible explanation is 
that the frequency of dying sources in cluster is higher Radio 

* The dying source B2 0924+30 is member of the nearby poor cluster 
WBL224 (White et al. 1999), but no extended X-ray emission is present 
in the RASS. 
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sources in dense environment, like e.g. the center of cooling core 
clusters, may have a peculiar accretion mode which results in a 
bursting duty cycle sequence of active and quiescent periods. Of 
the five dying galaxies presented in this work, WNB 1829+69 11 
and B2 0120+33 show evidence for a flat-spectrum core. A very 
faint core is detected also in the MERLIN full-resolution image 
of WNB 185 1+5707 a, while a steep spectrum feature is detected 
at the center of the host galaxy of WNB1734+6407. Thus, al- 
though the extended fading radio lobes in these sources are not 
powered by kpc-scale jets anymore, the AGN is still radio active 
at very low levels. 

To investigate these hypotheses we need to compare in detail 
the actual fading radio structures with the properties of the X-ray 
emitting gas. We recently have had the chance to observe three 
more of these clusters, Abell 2276, ZwCl 1829.3+6912, and 
RXC J 1 852. 1 +57 1 1 , with the Chandra satellite. We report on the 
results of these observations on a forthcoming paper (Murgia et 
al., in preparation). 

7. Summary 

In this work we presented the study of five 'dying' radio 
galaxies belonging to the Westerbork Northern Sky Survey 
minisurvey and to the B2 bright sample: WNB 1734+6407, 
WNB1829+6911, WNB1851+5707, B2 0120+33, and B2 
1610+29. 

These sources have been selected on the basis of their ex- 
tremely steep broad-band radio spectra, which is a strong indi- 
cation that these objects either belong to the rare class of dying 
radio galaxies or that we are observing 'fossil' radio plasma re- 
maining from a previous nuclear activity. We derived the relative 
duration of the dying phase from the fit of a synchrotron radia- 
tive model to the radio spectra of the sources. The modeling of 
the integrated spectra and the deep spectral index images ob- 
tained with the VLA confirmed that in these sources the central 
engine has ceased to be active for a significant fraction of their 
lifetime although their extended lobes have not yet completely 
faded away. In the cases of WNB 1829+69 1 1 and B2 0120+33, 
the fossil radio lobes are seen in conjunction with a currently ac- 
tive core. We found that WNB 1851+5707 is in reality composed 
by two distinct dying sources that originally appear blended to- 
gether at the lower angular resolution of the WENSS. 

All the dying sources of our sample are located (at least in 
projection) at the center of an X-ray emitting cluster 

Although no firm conclusions can be drawn due to the small 
number statistics involved, our results suggest that the duration 
of the dying phase for a radio source in clusters can be signifi- 
cantly higher with respect to that of a radio galaxy in the field. 
The simplest interpretation for the tendency of dying galaxies to 
be found in cluster is that the low-frequency radio emission from 
the fading radio lobes last longer if their expansion is somewhat 
reduced or even stopped by the pressure of a particularly dense 
intra-cluster medium. Another possibility is that the occurrence 
of dying sources is higher in galaxy clusters. Of the five dying 
galaxies, WNB1829+6911 and B2 0120+33 show evidence for 
a flat-spectrum core. A very faint core is detected also in the 
MERLIN full-resolution image of WNB 185 1+5707 a, while a 
steep spectrum feature is detected at the center of the host galaxy 
of WNB1734+6407. Thus, although the extended fading radio 
lobes in these sources are not powered by kpc-scale jets any- 
more, the AGN is still radio active at very low levels. Indeed, 
this may suggest that radio sources in dense environment, like 
e.g. the center of cooling core clusters, have a peculiar accretion 
mode which results in a bursting duty cycle sequence of active 



and quiescent periods. This result could have important impli- 
cations for theories of the life cycles of radio sources and AGN 
feedback in clusters of galaxies but awaits confirmation from fu- 
ture observations of larger, statistically significant, samples of 
objects. 

We may expect the existence of a large population of dying 
radio sources that have been missed from the current surveys 
because of very steep spectra. These sources are very faint at 
centimeter wavelengths but should still be visible at frequency 
below 100 MHz if they are only subject to radiative losses. Due 
to their sensitivity and angular resolution the upcoming low- 
frequency radio interferometers (such as LOFAR and LWA and, 
in the next future, the SKA) represent the ideal instruments to 
discover and study in detail these elusive objects. 
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Table 5. Integrated spectra. 
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r\ 1 c 1 r\ r\c\ 

0.35 ± 0.09 


1 .0 ± 0.3 


WNB 1851+5707 


38 


1800 ± 818 


0.818 


8C, Hales et al. 1995 


114± 11.9 


1686 ±818 




74 


1240 ± 150 


1.0 


VLSS, Cohen et al. 2007 


65.7 ±5.7 


1174 ± 150 




151 


723 ± 72 


0.904 


6CV, Hales et al. 1993 


36.3 ± 2.4 


687 ± 72 




232 


530 + 50 


1.0 


MIYUN, Zhang et al. 1997 


25.5 ± 1.4 


505 ± 50 




325 


387 ± 4.4 


1.0 


WENSS, RengeUnk et al. 1997 


19.3 ± 0.9 


368 + 5 




303 


A 1 A 1 
4j / ± 41 


1 f\A 1 

1.041 


ihiAAs, Douglas et al. lyyo 


C1 O 1 1 o 

51.0 ± 1.0 


OOC 1 /II 

3o5 ± 41 




408 


270 ± 10 


1.091 


"D 1 1 ..1 1 AAA 

d3.3, Altieri et al. 1999 


16.0 ± 0.7 


254 ± 10 




14UU 


3i ± l.i 


l.U 


JN Vaa, Condon et al. iWo 


C -70 1 A 1 O 

J. lis ± 0. 12 


/IT 1 O O 

4/ ± l.i 




1425 


41.9 ± 0.1 


1.0 


VLA, this work 




41.9 ±0.1 




4850 


6.1 ± 1.0 


1.0 


EfFelsberg 100-m, this work 


O AA 1 A AA 

2.U9 ± 0.09 


4.0 ±1.0 




4850 


O O C 1 A Al 

z.8j ± 0.03 


1.0 


VLA, this work 




O OC 1 A AO 

2.85 ± 0.03 




10450 


1.2 ± 0.3 


1.0 


EfFelsberg 100-m, this work 


1 1 O 1 A AO 

1.12 ± U.Oo 


< 1.0 


B2 0120+33 


74 


2847 ± 200 


1.0 


VLSS image, this work 


- 


2847 ± 200 




151 


1690 ± 170 


0.904 


6CV1, Hales et al 1993 


96.3 ± 3.6 


1594 ± 170 




235 


1100 ±88 


1.0 


GMRT, Giacintucci et al. in prep. 


- 


1100 ±88 




324 


727 ± 11.7 


1.0 


VLA, this work 


- 


727 ± 11.7 




325 


715 ± 16 


1.0 


WENSS image, this work 




715 ±16 




/I AO 

408 


£i£Lf\ 1 '^A 

660 ± 20 


1 A 

1.0 


heretti & Giovannmi 1980 


CA C 1 O C 

50.5 ± 2.5 


^CAA 1 OA 

609 ± 20 




610 


372 + 18.6 


1 A 

1.0 


GMRT, Giacintucci et al. in prep. 




OTO 1 AO £1 

372 ± 18.6 




1 Af\{\ 

1400 


1 AO 1 A C 

108 ± 0.5 


1 A 

1.0 


NVSS image, this work 




1 AO 1 A C 

108 ± 0.5 




1400 


99.5 ±1.5 


1.0 


VLA, this work 




99.5 ±1.5 




2380 


35 ± 3.5 


1.0 


AKbClBU, Uressel & Condon 19/8 




35.0 ± 3.5 




4860 


8.6 ±1.1 


1.0 


VLA, this work 




8.6 ±1.1 




8460 


< 3.1 


1.0 


VLA, this work 




< 3.1 


B2 1610+29 


74 


1055 ± 157 


1.0 


VLSS image, this work 


- 


1055 ± 157 




151 


774 ± 9.5 


1.237 


7C, Riley et al. 1989 


- 


774 ± 9.5 




325 


412 ±5.2 


1.0 


WENSS, Rengelink et al. 1997 


- 


412 ±5.2 




408 


349 ± 35 


1.091 


B2, CoUa etal 1970 


- 


349 ± 35 




610 


240 ± 12 


1.0 


GMRT, Giacintucci et al. 2007 




240 ± 12 




1400 


110± 11 


1.0 


NVSS image, this work 




110± 11 




1452 


101 ±2 


1.0 


VLA, this work 




101 ±2 




1502 


98.8 ±2 


1.0 


VLA, this work 




98.8 ±2 




2639 


49.8 ± 3 


1.0 


EfFelsberg 100-m, this work 




49.8 ± 3 




4710 


21.4 + 0.3 


1.0 


VLA, this work 




21.4 + 0.3 




4850 


19.0 + 0.4 


1.0 


EfFelsberg 100-m, this work 




19.0 ± 0.4 




8350 


4.8 ± 0.6 


1.0 


EfFelsberg 100-m, this work 




4.8 ± 0.6 



M. Murgia et al.: Dying Radio Galaxies in Clusters 
Table 6. Physical parameters of the dying sources. Uncertainties on fit parameters are at l-cr level. 



Source 




MHz 


foFp/ 


Ll51 

W/Hz 


^min 


^min 

erg/cm^ 


Myrs 


ta 
Myrs 


foFF 

Myrs 


WNR 1 734+6407 




230+1550 


3+"^ 

"•-'-0.2 


2.6 X 10^^ 


10 


8.3 X 10"'^ 


86 


60 


26 


WNB1829+6911 


"•'-0.2 


cn+1500 

-'"-20 




7.6 X 10^5 


7.8 


5.3 X 10-12 


218 






WNB1851+5707a 


"•-'-0.2 


220+620 
^^"-170 


0.2!°:? 


1.8 X 10^5 


9.8 


8.9 X 10-'2 


90 


72 


18 


B2 0120+33 Global spec. 
Sp. index profile 


0.6 (fixed) 
0.65 (fixed) 


440!?i° 
1840:t 


"•-'-0.1 

>0.8 


1x10^ 


5.2 


2.4 X 10-'2 


141 
69 


99 
< 13 


42 
>55 


B2 1610+29 Global spec. 
Sp. index profile 


"•"^-0.06 

0.6 (fixed) 


4770+f°« 


>0.39 

OKI 


1.4 X 10^ 


3.2 


9.3 X 10-" 


55 
58 


<33 
12 


>22 
46 



